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1. INTRODUCTION 

Massive multiple input multiple output (MIMO) is the key technology of fifth generation (5G) 
wireless communication systems. Theoretically, massive MIMO employs a large number of antennas up to 
thousands, but only several dozens or hundreds of antennas are practically implemented [1], [2]. In massive 
MIMO, the base station is generally equipped with a large number of antennas to serve one or multiple users 
where each user equipped at least with one antenna [3], [4]. Using massive MIMO can provide many benefits 
for 5G wireless systems by increasing the spatial multiplexing gain and diversity gain which lead to increase 
system capacity and link reliability, also enhance the spectral efficiency and energy efficiency [5], [6]. By 
increasing antennas number of massive MIMO system, the physical size of massive MIMO array will be 
large and practically unfeasible. To overcome this difficulty, massive MIMO system along with millimeter 
wave transmission (mmWave) are used. Smaller wavelengths enable operations with more antenna elements 
per system within small area size [6]-[8]. For example, the size of massive MIMO with (8*8) uniform 
rectangular array (URA) operates at a frequency of 6 GHz is approximately 300 cm?, while its size when 
operates at a frequency of 28 GHz is approximately 14 cm?. For mmWave massive MIMO system, the 


Journal homepage: http://beei.org 


1410 O ISSN: 2302-9285 


uniform planar arrays are more considerable so more antenna elements can pack in a compact-sized array and 
allowing beamforming in the azimuth and elevation domains (3D beamforming) [8], [9]. 

Operating at mmWave frequencies puts new challenges for 5G massive MIMO such as high path 
loss due to atmosphere attenuation, propagation challenges due to penetration, and RF chain complexity. The 
mmWave channels have low rank and it is specular. They tend to have a less multipath components 
compared to microwave, so the channel matrix is sparse with a smaller number of components due to the 
geometric nature of the mmWave channel [9], [10]. 

The propagation issue due to penetration can be contained through adopt small size 5G network cell 
such as microcell or femtocell. High path loss at mmWave band can be reduced by using digital array signal 
processing involved spatial multiplexing and beamforming. In spatial multiplexing, multiple data streams are 
encoded (e.g. alamouti coding) then transmitted independently using array antenna elements. At each 
receiver, the dedicated original signal can be reconstructed effectively [11]. 

Beamforming is a technique steers the transmitted signals of a massive MIMO system to the desired 
angular direction. Beamforming of massive MIMO array can be implemented either at baseband stage, which 
called digital beamforming or at radio frequency (RF) stage, which called analog beamforming [12]. The 
basic mechanism behind analog beamforming depends on using phase shifters to control the phase of each 
transmitted signal in order to steer the main beam into the desired angle [12], [13]. While in digital 
beamforming, the beamforming process is done using digital signal processing. Digital beamforming 
provides more flexibility and accuracy of the beamforming process in comparison to analog beamforming 
because it is easily to control the direction of beam through weight factors at baseband stage [14]. However, 
digital beamforming requires a separate transmit (receive) module dedicated to each antenna element, this 
module contains digital/analog converter analog to digital converter (ADC) or digital to analog converter 
(DAC), mixer, and power amplifier it is called RF chain. Unfortunately, design RF chain at mmWave is 
complex and expensive especially for massive MIMO array because a large number of RF chains is 
necessary [15]. 

Hybrid beamforming provides a trade-off between design complexity and performance. Hybrid 
beamforming is a combination of digital and analog beamforming [16], [17]. The concept of hybrid 
beamforming is first proposed in [18]. 

The essential process of hybrid beamforming is based on divided array antenna elements into many 
sub-array groups, only one RF chain component is attached to each group. Each array antenna element is 
connected to one phase shifter. This structure, significantly reduces the number of required RF chains, cost, 
and system complexity. The performance degradation of the hybrid beamforming can be controlled by using 
a well-designed adaptive beamformer. It is proved, if the number of RF chains is at least twice the number of 
data stream, the performance gap between hybrid beamforming, and fully digital beamforming approaches 
zero [19], [20]. 

The analog precoder maps the outputs of RF chain to the antenna elements. According to this 
mapping, there are two main structures for hybrid beamforming which are fully-connected structure and 
partially-connected structure. In fully-connected structure all array antenna elements are connected to each 
RF chain, while in partially-connected structure each RF chain is connected only to a set of the antenna 
elements [20]-[22]. Note that the term "precoder" means beamformer at the transmitted side, while the term 
"combiner" means beamformer at the receiver side. 

While the digital beamformer modify both the amplitude and phase of transmitted symbols, the 
analog beamformer has a constant amplitude 1/v (N,), and N-resolution quantized phases in fact. The 
resolution of analog phase shifters is limited by the number of bits used in it. So, high resolution phase shifter 
provides fine resolution analog beamforming but increases power consumption and hardware complexity 
[23], [24]. 

In this paper, a hybrid precoder design is proposed for a MIMO orthogonal frequency division 
multiplexing (MIMO OFDM) system based on the second order channel statistics (transmit covariance 
matrix) to reduce the overhead information. The joint spatial division and multiplexing (JSDM) scheme is 
modified to design a precoder serves users with multiple antennas. Avoiding the coupling between the analog 
precoder design and the analog combiner design (the combiner design is beyond the scope of this paper) is 
achieved according to the approach of Kronecker channel model [22], [25]. A wideband precoder for 
wideband mmWave channel is obtained by modelling this wideband channel as multiple narrowband 
subcarrier channels, so the proposed narrowband precoder is applied to each OFDM subcarrier channel. 
Although JSDM does not provide an optimal beamforming, but spectral analysis shows a little degradation 
between the proposed precoder and fully digital one especially when large number of antennas used at the 
base station. 

Notations: in this paper, the following notations is used, 
(OT, (OEC I Ile Il liz and EL ] denote the transpose, conjugate transpose, conjugate, frobenius norm, 
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12-norm, and statistical expectation respectively. CN (u, K) is a circularly symmetric complex Gaussian 
distribution with mean u and covariance matrix K. X®Y is the Kronecker product of X,Y matrices. 
Span(V), Spant(V) denote the column space of (V), and its orthogonal complement respectively. 


2. SYSTEM MODEL 

Consider a downlink multi-user multi-stream mmWave massive MIMO system operates at 28 GHz, 
where the base station (BS) equipped with a uniform planar array (UPA) of N, antennas, and communicates 
Ns independent data streams with Ny users each equipped with N,,; antennas, where i = 1, 2, ...... > Nu. Nsi 
represents the number of independent data streams per user, Np is the number of RF chain at the base 
station, where Nç = pA N,;, and N, = ya N, i- The users distributed within a range of 1000 meters from 
the base station. working at the mmWave frequency provide a narrower beam width, since the beam width is 
inversely proportional to the operating frequency, and this can enhance the beamforming process. Also, the 
large available mmWave bandwidth increases the transmission bit rate [8], [22]. However, a wideband 
beamforming is needed to take advantage of large mmWave channel bandwidth. OFDM technique is used to 
divided the wideband mmWave channel into many narrowband channels, which help in the design of 
beamforming process [26]. The block diagram of the system is shown in Figure 1. 
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Figure 1. Massive MIMO-OFDM system block diagram 


Hybrid beamforming is applied at the base station, the digital beamformer is used for each OFDM 
subcarrier and the analog beamformer is common to all OFDM subcarriers. The hybrid beamforming weights 


t 
of the k* OFDM subcarrier is of form TE = Ferri fei where Fpr i E C™E*NRF.i is the analog precoder and it 


is common for all the subcarrier, and fini e CNRriNsi is the baseband precoder of the k™ subcarrier. In 
spite of that the users are equipped with multiple antennas, but no hybrid combiner is applied here, so the 
paper is focused only on the design of the hybrid precoder based on a modified JSDM for a wideband 
mmWave MIMO-OFDM system. However, it has been shown later that the design of hybrid precoder and 
combiner are independent of each other, also the hybrid combiner design can use the same model proposed 
here. 

For each OFDM subcarrier, the symbol signal (xk € Cs.i*!) which is belong to 7” user is precoded 
using digital baseband precoder fist then the precoder output is passed through OFDM, then is upconverted 
to RF stage using Néri RF chains. The fully-connected structure is used to connect the output of RF chains to 
the elements of the base station antenna array, so the analog precoder Fpp, uses Neri x N, phase shifters. 
According to the above process, the transmitted signal (y¥ € C1) is: 


Yi = FeeifepiXt (1) 
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The received signal (rë € Cris) at i user is: 


= JP; HE FrrifgeiXt +e VP; HE Frr,j fle jX +N; (2) 


j#i 


P; is the received average power at the i® user, n; E C%ri*? is the noise signal. n; has complex gaussian 
distribution with zero mean and covariance matrix Rn; = E[n;n'7] (e.g., ni ~ CN (0, Rn;)). HE E CNritNe 
is the channel matrix between the base station and the i user at k" subcarrier. The Saleh- Valenzuela channel 
model [27] which adopted in this work to modelling the propagation of MIMO channel is: 


Nei 


Nee, 
Sy Xe F D 1 %emitri(Ppmir Omi) af (Pimi bmi) (3) 


This narrowband channel model assumed to be a sum of N,g,; rays divided into Nop; clusters, and each 
cluster contains Np; rays for £ = 1 ,2..., Negi Y = J NeNr/Nrayi> Nray,i = y Nei- emi represents the 
complex gain of the m" ray in the #*” cluster and it is assumed to be zero mean independent random variable. 
mi 92m, are the azimuth and elevation angles of arrival (AoA) respectively, which are random variables 
with mean Pri OF i. Pimi Ob mi are the azimuth and elevation angles of departure (AoD) respectively, which 
are random variables with mean Pei 96. The AoD and AoA are assumed to be independent random 
variables. Finally, the a; (bimi Femi) and a, ; (bimi emi) are the antenna array response vectors of the 


base station and the i" user respectively. The channel model of the above equation is modified for wideband 
MIMO-OFDM mmWave channel model as in [28], so 


a Nei 


HE 


Èm Bemkiâri (bimi O3 mijat (Dimi Oo mi) (4) 


The right hand of the above equation is similar to the narrowband channel model of (3), where Bemxi = 
Aem ipek i represents the complex gain of the m™ ray in the € cluster for k*” subcarrier, k = 1,2,..., Ngc- 
Peri = E? pnt, — Tp) eJ2mkn/Nsc, Nyc is the number of OFDM subcarrier, and Nop is the number of 
cyclic prefix (CP) and should be greater than the largest time delay of the channel to remove the intersymbol 
interference (ISI). y(t) represents the function of transmitter (receiver) lumped pulse shaping. T, is the sampling 
time, and Zp is the average time delay of the ” cluster. vec (HE)~CN (0, KF) and KE E CNNr NN ri 
represent the channel covariance matrix. Hence, each subchannel H¥ can be considered as a narrowband Saleh- 
Valenzuela channel model, so any narrowband beamformer design algorithms can be implemented for each 
subchannel HE. Moreover, the Kronecker channel model [Ref. 170, Ref. 84] assumes that the transmit and 
receive correlation are separable, based on that assumption, the channel covariance matrix K¥ can be expressed: 


KF = KE Q KA (5) 


H H f : : 3 
Where, Kf; =E [ack nk] ,and Kh, =E | sekok | are the transmit and receive covariance matrices of 


the channel. based on (5), the design of analog precoders and combiners can be independent of each other. 
Without loss of generality, the following assumption are made: equal power distribution over all data 


SPEDEN x i : H 1 
streams, x¥ has a normal distribution with zero mean covariance matrix Rx; = E [acct | = go Nsi (e.g., 
si , 
x; ~ N (0,Rx;)). rË transmitted through a narrowband block fading channel. This assumption can be 
feasible by using OFDM technique which converts wideband selective fading channel into several 


narrowband flat fading channels [28], [29]. {||Ferifse.i j = Nsí V k, i} to satisfy the total power constraint. 


3. PROBLEM FORMULATION 

The main step of the hybrid beamforming design is to find the optimal values of precoder weights 
( fae Frei) that minimize the MSE between xE and rë over all users and subcarriers. The problem can be 
written as: 


min ONY 8050 [ak — rë l] 
pFRF,i 
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s.t. Peil = Nov Vi © 
Frei © Fre 


A popular hybrid beamforming design method is JSDM is considered in this work. JSDM is 
designed for multiuser, each one with single antenna. JSDM is used to provide beamforming coefficients at 
base station only [30]. 

In this work, a modified version of JSDM is used to obtain precoder weights for multiuser, each 
with multiple antennas. For simplicity, one user per group is assumed and the users are well spatial separated. 
First, the analog precoder is designed based on the second order channel statistics (channel covariance 
matrix) and not on the instantaneous channel state information which reduce the overhead information needs 
to process and fed back because channel covariance matrix is varying slowly with respect to coherence time. 
At the beginning, each user estimates its channel and sending V# to the BS. the singular decomposition value 
(SVD) of channel matrix (HE) is defined as: 


dk = yey (7) 


Where UF € CNni*Nri, V¥ e CNt*Nt are the left-singular vectors and right-singular vectors respectively. 
xe e RNni*Ne is a diagonal matrix with rank (Zf) = N,;, and only non-negative real elements and these 


f : F : H 
elements represent the square root of the eigenvalues of transmit covariance matrix KË = E [74K HI. 


Furthermore, the right-singular vectors matrix V;* can be obtain from eigenvalue decomposition (EVD) of 
transmit covariance matrix is: 


2. KH 
KE = VEE Vi (8) 
Following the approach of [28], an approximate of EVD of K¥ is defined as: 
Ki = y’V AFV? (9) 


Where Af € C^ri*Nri is a diagonal matrix, whose elements are o2,,\Pexil Ae. Ay, is a diagonal matrix 
contains the dominant eigenvalues of transmit covariance matrix, and V; E€ CNt*Nri is a semi-unitary matrix 
contains the eigenvectors corresponding to the dominant eigenvalues of Ay. Note that only A¥ depends on k* 
subcarrier and V; does not. V; is an approximate of the first N, į eigenvector columns of VE Based on the 
approach of [30], [31], =; is defined at the base station as: 


_ T 
Si = [P V2 ie Vi-1 Vint ooo Vuy-1 Vya] y 


Where V, is a complex conjugate of the first N,,, columns average of V** over all subcarriers, 2; matrix 


contains the eigenvectors of all users except that of i user. So E; € CEN: and = ai N, j. The 
SVD of &; is: 

E; = WEB” (11) 
B; E CNe*Ne and its rank=L;. B; can be written as: 

B= [g g] (12) 


Where, Bi € C%*4i contains the right eigenvectors corresponding to the dominant eigenvalues of E;, and 
BO e CNex(t-£i) contains the right eigenvectors corresponding the weakest eigenvalues of =;, which 
represents the null space of =; such that Span (8; 0) = Spant d3: Vi + i}). The effective channel Hk € 
CNrix(Ne-£i) is defined as: 


FE = HEBO (13) 
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The analog precoder Fp; ; can be obtained by concatenating the right eigenvectors corresponding to the 
dominant eigenvalues of the effective channel H* with the projection onto Span(B). Ŭi; is the columns 
average of the eigenvectors V* over all subcarriers, where V;* is obtained from the eigenvalue decomposition 
of the transmit covariance matrix KE = ese ye Rewrite P, C™ t-LDXN Ri as: 


K=O 9] (14) 


; 5 (1) ih 
As mention above, the product of Band Vi — produces Neri column vectors that maximize the 
transmission rate and eliminating inter-group interference for i user [32]. These column vectors represent 
the analog precoder Fp, ;. 


= (1) 
Feri = 8; * D, (15) 


Based on (15), the second term of the right-hand side of (2) can be eliminated since the effective 
channel (É 5 = HEBO = 0) for exact block diagonalization (BD), so (2) becomes: 


rë = JP; Hi Fer iff xk + + ni (16) 


Moreover, span (V;) £ span ({V; :j # i}), dim(span(V,) N span ({V; :j # i}) = Ns; should be satisfied 
to achieve exact BD, and the condition span (Fri) & span* dv JF i}) is necessary for zero-forcing. The 


design of hybrid precoder is summarized as shown in. Given Fp, ;, the subcarrier digital precoder fav can be 
xN 


= Nri 
obtained by define the effective channel AE € C NRF as: 


RE = HF Far; (17) 


Algorithm 1 Analog Precoder Fpr; Design 


1- Estimate H* at each user. 

2- Define KK = E[H¥" HF]. 

3- Compute the EVD of (Kj) to find the eigenvectors (V;*). 
4- Each user feedback its eigenvectors (VX) to the BS. 
5- Define & according to (10). 

6- Find gO according to (11), and (12). 

7- The BS broadcasting B,O for users. 

8- Define the effective channel CHS according to (13). 
9- Compute the EVD of (KE) to find the eigenvectors (Vj). 
10- Find DB,” according to (14). 

11- Finally, the analog precoder is defined according to (15). 


Based on that Far’ Frei, = Clyt for large BS antennas N, [7], [28], the subcarrier digital precoder Jepi is 


defined as the norm-normalized N,; dominant eigenvectors of covariance matrix Fp ve KK Fp Fiy So, the 
subcarrier precoder is the product of analog RF precoder by subcarrier baseband precoder. 


4. SIMULATION RESULTS 

An evaluation of the proposed system via simulation is presented in this section. A multi-users 
multi-streams OFDM massive MIMO system, shown in Figure 1, is considered. The users are randomly 
distributed within 1000 m diameter area around the base station in a manner ensure a well spatial separation 
between users. The base station communicates with the users over multi-clustered channel as in (4) with ten 
clusters each with ten rays, so the total number of rays is one hundred. The system using 256 subcarriers 
OFDM with 64 cyclic prefix to handle the wideband mm wave (28 GHz) channel. Channel sounding is used 
to estimate the channel by sending a preamble signal. Without loss of generality, each user is equipped with 
eight antennas and two data streams is dedicated for each user to ease the calculations. 

First, the result of the proposed system is compared to that of fully digital precoding system with 
full channel state information for a variety number of base station antennas N, = 64,128,256 and 512. The 
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antennas are arranged in UPA. Figure 2(a) shows the average sum rate versus signal to noise ratio (SNR) for 
the proposed and fully digital system. The number of users is eight, so the total number of receivers’ antennas 
is 64 therefore the minimum number of base station antennas must be more than 64. As shown in Figure 2(a), 
the gap between the proposed system and fully digital precoder system reduce as the number of transmitted 
antennas is increased, while the gap at O dB is about 8 for 64 antennas it does not exceed 0.7 for 512 
antennas. The large degradation in the proposed system when the transmit array antennas is 64 is due to that 
the number of transmit array antennas should be more the total number of array antennas of all receivers, 
which is unverified condition here. Therefore, when the number of users is six and the total receive array 
antennas become 48, the gap for 64 transmit antennas is enhanced significantly in comparison to that for 
higher transmit antennas as shown in Figure 2(b). In Figure 3, the system is tested for the quantization effect 
of the phase shifters of the analog precoders. When using 1-bit phase shifters instead of unquantized phase 
shifters, the spectral efficiency degradation for 512 transmit array antennas is about 1.55 at O dB while it is 
about 0.43 for 128 transmit array antennas. 
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Figure 2. Spectral efficiency of the proposed system for (a) 8 users (b) 6 users. Each user equipped with 8 
antennas, 2 data stream per user, and number of RF chains is 48 (three chains per data stream) 
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Figure 3. Spectral efficiency of the proposed system for 8 users each equipped with 8 antennas, 2 data stream 
per user, and number of RF chains is 48 (three chains per data stream), considering quantization effect on 
analog phase shifters (1 bit phase shifters) 


A comparison of average sum rates at O dB for different number of transmit antennas when the 
number of users is 4 and 8 are listed in Table 1 considering the quantization effect of phaser shifters. The 
quantization of phase shifters has a limited effect on the system because the JSDM algorithm is not an 
optimal solution for the design of hybrid precoder. Figure 4 shows the effect of number of users on the 
average sum rates for different numbers of transmit antennas at 0 dB signal to noise ratio. 
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Table 1. Average sum rates for different transmit antenna number at 0 dB 
Unquantized phase shifter 1-bit phase shifter 2-bit phase shifter 3-bit phase shifter 


Nt=64 7.661 7.409 7.566 7.656 
i Nt=128 11.578 10.850 11.395 11.535 
Z Nt=256 16.091 14.796 15.539 15.856 
Nt=512 20.452 18.735 19.721 20.126 
co  Nt=128 8.791 8.357 8.699 8.732 
i Nt=256 13.147 12.043 12.704 12.974 
Z Nt=512 17.660 16.102 16.997 17.350 
25 f —x— N =64 
—o— N5128 


Average sum rate (bps/Hz) 


Number of Users 


Figure 4. Average sum rates versus number of users for different transmit array size at (SNR=0 dB) 


5. CONCLUSION 

In this paper, a hybrid precoder is designed for multi-user multi-stream wideband MIMO-OFDM 
system based on JSDM for downlink transmission. The effect of transmit array size on the proposed system 
is discussed and the condition that the antennas number of transmit array must be more than the antennas 
number of all users array is verified. Also studied how sum rates effect by limiting the number of bits of the 
analog phase shifters which has small effect on the system because JSDM do not provide an optimal 
precoder, so the effects of the number of the phase shifter bits is limited. Even though the JSDM does not 
provide an optimal precoder but the results show a little degradation in the spectral efficiency. The proposed 
system provides a good computational efficiency because it does not depend on the instantons channel state 
information, also this will reduce the overhead information and the time to fed it back. 
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